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Abstract 

We study the Hawking process on lattices falling into static black 
holes. The motivation is to understand how the outgoing modes and 
Hawking radiation can arise in a setting with a strict short distance 
cutoff in the free-fall frame. We employ two-dimensional free scalar 
field theory. For a falling lattice with a discrete time-translation sym- 
metry we use analytical methods to establish that, for Killing fre- 
quency lo and surface gravity k satisfying k <C w 1 / 3 < 1 in lattice 
units, the continuum Hawking spectrum is recovered. The low fre- 
quency outgoing modes arise from exotic ingoing modes with large 
proper wavevectors that "refract" off the horizon. In this model with 
time translation symmetry the proper lattice spacing goes to zero at 
spatial infinity. We also consider instead falling lattices whose proper 
lattice spacing is constant at infinity and therefore grows with time 
at any finite radius. This violation of time translation symmetry is 
visible only at wavelengths comparable to the lattice spacing, and it 
is responsible for transmuting ingoing high Killing frequency modes 
into low frequency outgoing modes. 
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1 Introduction 



As field modes emerge from the vicinity of the horizon they are infinitely 
redshifted. In ordinary field theory there is an infinite density of states at 
the horizon to supply the outgoing modes. How do these outgoing modes 
arise if the short distance physics supports no infinite density of states? And 
how does the short distance physics affect the Hawking radiation in these 
modes? By insisting on a fully sensible resolution of the apparent conflict 
between black holes and short distance finiteness we hope that some deep 
lessons can be learned. That is the underlying motivation for the present 
work. 

One way to avoid an infinite density of states is to have some physical 
cutoff at short distances, related to a fundamental graininess of spacetime.0 
Analogies with condensed matter systems such as fluids, crystals, Fermi liq- 
uids, etc. suggest that in this case the long wavelength collective modes 
which are described by field theory will propagate with a dispersion relation 
that deviates from the linear, Lorentz-invariant form at high frequencies. A 
number of (two dimensional) linear field theory models with such behavior 
have now been studiedjl], |2], ^, ||. It turns out that the Hawking radiation 
is extremely insensitive to the short distance physics, as long as neither the 
black hole temperature nor the frequency at which the spectrum is examined 
is too close to the scale of the new physics.0 What is striking is that this is so 
even though the behavior of the field modes is rather bizarre: the outgoing 
modes that carry the Hawking radiation arise from exotic ingoing modes that 
bounce off the horizon.f] 

3 String theory provides a different way, in which the states are iess localizabie than in 
ordinary field theory. 

4 The ultra low frequency part of the spectrum, which has not yet been computed, 
might turn out to be non-thermal. 

This is what happens in the case where the group velocity is subluminal at high 
frequencies. In the superluminal case the outgoing modes arise from superluminal modes 
that emerge from behind the horizon. Ultimately these modes come from the singularity 
(for a neutral black hole), so it is not so clear one can make sense of this case. However, 
UnruhJ^] and CorleyQ have recently shown that if one simply imposes a vacuum boundary 
condition on these modes behind the horizon the usual Hawking radiation outside the black 
hole is recovered. A black hole with an inner horizon (such as a charged one) behaves very 
differently in the superluminal case however. The ergoregion inside the black hole is 
unstable to self-amplifying Hawking radiation[Q. 
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Although these models do provide a mechanism for generating the out- 
going modes without an infinite density of states at the horizon, they still 
behave unphysically: wavepackets cannot be propagated backwards in time 
all the way out to "infinity" (i.e. the asymptotic region far from the black 
hole). For example, using Unruh's dispersion relation]]]], which has a group 
velocity that drops monotonically to zero at infinite wavevector, the wavevec- 
tor diverges as the wavepacket goes (backwards in time) farther from the 
black hole. So, in this case, the infinite redshift is just moved to a new lo- 
cation. Evidently the theory is pushed into the trans-Planckian regime after 
all. To make sense of — or at least to sensibly model — the true origin of the 
outgoing modes, it therefore seems necessary to work with a theory that has 
a well-behaved physical cutoff. A simple way to implement such a cutoff is 
to discretize space, and work with a lattice theory preserving continuity in 
time. This is roughly similar to what is happening in a condensed matter 
system, but we can preserve strict linearity for the lattice field theory and 
still model the key effect of the cutoff. 

In this paper we study two lattice models of this nature obtained by 
discretizing the spatial coordinate in a freely falling coordinate system. (If 
instead a static coordinate is discretized then the lattice points have diverging 
acceleration as the horizon is approached, and their worldlines are spacelike 
inside the horizon. This leads to pathological behavior of the field.) If there 
really is a cutoff in some preferred frame, then that frame should presumably 
fall from the "cosmic" rest frame at infinity in towards a black hole. This 
would be like Unruh's sonic analog of a black hole || [T], [| or the helium- 
3 texture analog pTOfl , where the short distance cutoff is provided by the 
atomic structure of the fluid which is freely flowing across the phonon or 
other quasiparticle horizon. 

One particular choice of discretization has the feature that a discrete 
remnant of time translation invariance survives. This makes the model easier 
to study analytically, and we exploit this to show here in Sections that 
in a leading approximation the black hole radiates thermally at the Hawking 
temperature. The same result was found previously by Unruh |TTJ using 
numerical evolution of the lattice field equation. Unfortunately, however, 
this particular lattice model is still not satisfactory as a model of physics 
with a fundamental cutoff, because the proper lattice spacing goes to zero at 
infinity! 

It is easy to avoid the vanishing lattice spacing by discretizing instead 
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a spatial coordinate which measures proper length on some initial spacelike 
surface all the way out to infinity. However, since we also want the lattice 
points to fall freely into the black hole, this results in a lattice spacing that 
grows in time, as shown in Section |6|. The growth of the lattice spacing 
suggests that we have still not found a satisfactory model with a short dis- 
tance cutoff. (An alternative which avoids this problem will be discussed at 
the end of this paper.) However, it is rather instructive to understand the 
physics of this model with the growing lattice spacing. In this model time 
translation symmetry is violated for short wavelength modes but not for long 
wavelengths. In fact, the Killing energy of an outgoing mode can be much 
lower than that of the ingoing mode that gave rise to it. This is essential 
to producing the outgoing long wavelength modes in this model since a long 
wavelength ingoing mode will of course sail across the horizon into the black 
hole rather than converting to an outgoing mode. This mechanism is studied 
in Section ^| with the help of the eikonal approximation. 

We adopt units in which h — c — 5 — 1, where 5 is the coordinate lattice 
spacing, and we use the "timelike" metric signature. 



Our goal is now to "latticize" the theory of a scalar field propagating in a 
static black hole spacetime. For each spherical harmonic, the physics reduces 
to a two-dimensional problem in the time-radius subspace. The short dis- 
tance phenomena we wish to study have nothing to do with the scattering 
of modes off of the angular momentum barrier, so nothing essential is lost in 
dropping the angular dependence and studying instead the physics in a two 
dimensional black hole spacetime. 

We begin with a generic static two dimensional spacetime, and choose 
coordinates so that the line element takes the form 



The curves of constant x are orbits of the Killing field \ — &t- The curves with 
dx = v(x)dt are geodesies which are at rest with respect to the Killing field 
where v (x) = 0, and the proper time along these geodesies is t. The constant 
t time-slices are orthogonal to these geodesies, and the proper distance along 



2 Falling lattice models 




(i) 
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these time-slices is x. In appendix [A] we explain why such a coordinate system 
can always be chosen. 

To represent a black hole spacetime with an asymptotically flat region 
at x — > oo, we choose v{x) to be a negative, monotonically increasing, 
function with v{oo) = 0. The event horizon is located where the Killing 
vector becomes lightlike, i.e. where v(x) = — 1. For the Schwarzschild 
black hole this coordinate system corresponds to the Painleve-Gullstrand 
coordinates [fT2|, |13|| , with x = r and v(x) = —J2GM/x. A sketch of the re- 
lation between these coordinates and the ingoing Eddington-Finkelstein null 
coordinate v is given in figure |l[ (The wavepacket trajectory is discussed in 
section II.) 



horizon 



t = 

constant 



ingoing 
ughtray 



wavepacket 




Figure 1: Painleve-Gullstrand coordinates and ingoing light rays. The trajectory 
of a wavepacket that is outgoing with low wavevector at late times is sketched. 

A new coordinate y that is constant on the free-fall worldlines dx = v(x)dt 
is defined by 

» : = '-/!)■ (2) 

This yields the line element 

ds 2 = dt 2 - v 2 (x)dy 2 (3) 
where x is now a function of t — y obtained by solving (0) for x. In these 
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coordinates the Killing vector is given by 



X = dt + d y . (4) 
The action for a real scalar field in these coordinates is 

S = l -Jdtdy^g^d^d^ (5) 

dtdy(\v{x)\{d t 4>Y- —L^) 2 ^ (6) 

and the equation of motion is 

-*+ + <^-«'W + $S)) a <* = 0. (7) 

We could now alter the theory to include high frequency dispersion by 
replacing d y by F(d y ) = d y + ady + ... in the action (|). This is similar to what 
was done in the models already studied]!], |2], |3|, the only difference being that 
there it was d x that was replaced by F(d x ). Since d y = —v(x)d x , these two 
modifications are essentially the same near the horizon where v(x) — — 1, and 
in fact they are quite similar in all regions where v(x) is of order unity. It is 
only asymptotically, where v(x) goes to zero, that their behavior should differ 
substantially. We previously preferred to modify d x since it is the derivative 
with respect to proper distance on a constant t surface everywhere. Now 
however we want to discretize the spatial coordinate and, as explained in the 
introduction, we do not want to discretize x because it is infinitely accelerated 
at the horizon. Instead, we discretize the free-fall coordinate y. 

One possible spatial discretization of the action is0 



where D is the forward differencing operator D(fi m (t) := (4> m+ i(t) — <fi m (t))/5, 
5 is the lattice spacing in the y coordinate, and v m (t) := v (x(t — m8)). In 
the remainder of this paper we shall work in units of the lattice coordinate 

6 For later convenience in the WKB approximation we take the average of v m+ \ and 
Vm in the second term in the action. 
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spacing, so that 5 = 1. Varying the action (^) gives the equation of motion 
for (p m (t), 

This lattice action has a discrete symmetry 

{t,m) -> (t + l,m+l) (10) 

which is the remnant of the Killing symmetry generated by (f|) . The meaning 
of this is that shifting forwards in time by one unit at fixed static coordinate 
x is just enough time for the next lattice point to fall from x(t, y + 1) to 
x(t,y). This symmetry will be heavily exploited in the following analysis.[] 

Note that the y coordinate is infinitely bunched up as v — > (see (|3|)), 
which occurs at infinity for a black hole type metric. Therefore the uniform 
discretization y m = m yields a proper lattice spacing that goes to zero at 
infinity. This is undesirable from a physical point of view, but it is a con- 
venient choice mathematically, since unlike other discretizations it preserves 
the symmetry fllPl) . Also, as long as we do not try to evolve the scalar field 
modes all the way to infinity, the decreasing proper lattice spacing is benign 
and has no effect on the physics of the Hawking process. However, since 
our goal is to understand how the outgoing modes can be accounted for in a 
theory that has a "reasonable" short distance cutoff, we shall return to this 
issue in Section ||. 



The lattice model defined by (Q) was studied numerically by Unruh[|TI 
He found by propagating wavepackets backward in time that the outgoing 
modes come from exotic ingoing modes and, if these ingoing modes are in 
their ground states, then the outgoing modes are thermally occupied at the 
Hawking temperature. In the next three sections we use analytic methods 
to understand the propagation of these wavepackets and the computation of 
the flux of radiation from the black hole. Our results are in agreement with 
Unruh's numerical results. 



7 The existence of this discrete remnant of the Killing symmetry was pointed out to us 
by W.G. Unruh. In Sections ^| and ^ we study a similar model in which a reparametrization 
of y is discretized and no discrete symmetry survives. 
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3 Lattice dispersion relation 



Due to the symmetry fllPf) of the lattice action (§) there exist mode solutions 
of the form 

(t> m {t) = e-^f{m-t). (11) 

Under the discrete symmetry (|i~0D the mode (|TT]) changes by a phase factor 
as <p m (t) — > e~ lu} 4> m (t) . This identifies uj as the Killing frequency which is 
defined modulo 2im and is conserved. 

To derive the dispersion relation we plug the ansatz 

(f) m {t) = e -i^ e iKm-t) = e -i(u>+k)t e ik m ^ 

into the equation of motion and treat v m (t) as a constant. The result is 

\v\(u + k) = ±2sin(Jfe/2). (13) 

The free-fall frequency, i.e. the frequency measured along the free-fall lines 
of constant y, is defined by dt<p = —iu$(f). The form of the modes fll2|) then 
shows that 

uj s = uj + k. (14) 

To understand what range of uj and k are considered distinct, note that 
the modes defined by flTJD are invariant under the simultaneous shifts 

k^k + 2nn (15) 
uj ^ uj — 2nir (16) 

for any integer n. Thus we can transform any (uj, k) pair into an equivalent 
pair (uj', k') where k! lies within a fixed range of length 2i\ (the standard 
choice being —it < k' < it). The value of uj' is unconstrained with this range 
of k! . One choice of fundamental domain of (uj, k) pairs is therefore given by 

— 7T < k < IT, — OO < UJ < oo. (17) 

Conversely, we could just as well use the above transformation to force uj' to 
lie within a fixed range of length 2tt leaving k' arbitrary. 

The dispersion relation (p~3[ ) has a useful graphical representation (see 
Fig. H): On a graph with abscissa k, the straight line with slope \v\ and 
fc-intercept — uj intersects the curve ±2sin(/c/2) at a A; that is a solution or 
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Figure 2: Graphical representation of the dispersion relation (13). 



"root" of the dispersion relation. A wavepacket constructed from modes 
of the form ([11]) with Killing frequency near u will propagate through the 
lattice spacetime with conserved Killing frequency. This propagation can be 
represented graphically in the WKB approximation by following a point on 
the dispersion curve. Since the Killing frequency is conserved, the fc-intercept 
of the straight line is fixed, while the slope \v(x) \ of the straight line changes 
according to where the wavepacket is located. The direction of motion with 
respect to the static position coordinate £ := y — t = m — £ is determined by 
the group velocity d£/dt which is given by 

v g = dw/dk = ± C ° 8 ^J 2) -1. (18) 

Therefore the sign of the group velocity is the sign of the difference between 
the slope of the ±sin curve at the intersection point and the slope of the 
straight line. The group velocity in terms of y is dy/dt = ±cos(fc/2)/|i>|, 
which is always less than the speed of light according to the line element (131). 
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4 Origin of the outgoing modes 



In this section we argue using the dispersion relation that outgoing low 
wavevector wavepackets indeed originate as ingoing high wavevector wavepack- 
ets which "bounce" off of the horizon. A spacetime diagram of the process 
is sketched in Figure [I]. 

To see where the outgoing modes come from, consider a late-time, positive 
Killing frequency, outgoing packet centered on a small positive wavevector 
k +s . This wavepacket is represented on the dispersion curve in Fig. [2] as the 
point labeled k +s . Following this back in time using the graphical method 
described in the previous section we find that it moves up the dispersion 
curve until it reaches the tangency point k tp at which the group velocity (|18|) 
vanishes. This is the turning point, where the WKB approximation fails. If 
uj <C 1, the straight line is extremely close to the sine curve for many k val- 
ues. This means that when the wavepacket is close to the horizon it is really 
a superposition of many k values, including negative ones. The amplitude of 
the negative wavevector piece, which determines the Hawking radiation, is of 
order exp(— ttuj/k) where k is the surface gravity of the horizon. The positive 
and negative wavevector pieces both propagate back away from the horizon, 
evolving into the modes k + and k_ respectively.^ Thus we see that the out- 
going positive Killing frequency modes come from ingoing large wavevector 
modes which "bounce" off the horizon. This continuous evolution from one 
type of modes to another is called mode conversion. The same phenomenon 
occurs in the continuum models in which the high frequency dispersion is 
put into the theory by adding higher spatial derivative terms to the action. 

Now let us compute the values of the wavevectors k± and k +s correspond- 
ing to a fixed frequency u as \v | — > at infinity. From the dispersion relation 
( |13D or Fig. H, one sees that all three wavevectors k +s , k_, and k + converge 
to zero modulo 2tm independent of the value of u. This rather strange result 

8 Other modes get excited as well, but only slightly. From "reflecting" off the back- 
ground curvature a small negative wavevector piece will arise. This will have extremely 
small amplitude however, for the following reason. There is no scattering at all for a mass- 
less scalar field in the continuum due to conformal invariance of the action. On the lattice 
this symmetry will remain approximately for wavelengths much longer than the lattice 
spacing, and short wavelength modes will not see the curvature. As v(t — m) becomes 
smaller, there are also more wavevector roots to the dispersion relation with |fc| > 2tt 
which are also presumably excited slightly by scattering. 
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follows because the continuum metric has the form ds 2 = dt 2 — v 2 dy 2 , and 
so the ^-lattice spacing goes to zero as v goes to zero. Therefore any mode 
of finite proper wavelength will have infinite coordinate wavelength and zero 
coordinate wavevector. To resolve these modes we can look at their proper 
wavevectors k p = k/\v\ instead of the coordinate wavevectors. For the k +s 
wavevector, as v — > we may approximate 2sin(A;/2) « k in the disper- 
sion relation fll3|) (with the plus sign), which yields k ~ \v\u, so the proper 
wavevector goes to just to. For the k± wavevectors, we first use the symmetry 
relation (|IqJ16l) to shift the coordinate wavevectors (and therefore also the 
frequency uS) so that they converge to zero as v — > 0, and then use the small 
k approximation in the dispersion relation to obtain the proper wavevectors 
k Pt ± = — (a;±27r). Therefore the late time, long wavelength, outgoing Hawk- 
ing particle arises from a pair of short but finite proper wavelength ingoing 
modes. It follows from the discussion below (F3) that, at spatial infinity, the 



group velocity for these wavevectors is equal to the speed of light. 

In the next section we compute the amplitudes of the k + and fc_ pieces of 
the ingoing wavepacket. Crucial to the validity of the approximation used in 
this calculation is the maximum value of the wavevectors in the wavepacket 
solution near the horizon. We can estimate this maximum by a simple cal- 
culation using the dispersion relation. The classical turning point is located 
where the straight line of figure |2] is tangent to the sine curve, labeled k tp 
in the figure. Although the wavepacket tunnels beyond the classical turning 
point, it is not propagating there, so its shortest wavelength near the horizon 
should be roughly given by the wavelength at the classical turning point. 
The wavevector at this point satisfies the dispersion relation (|13|) (with the 
plus sign) and the relation 

\v\ = cos(A;/2) (19) 

expressing equality of the slopes of the two curves. If uj <C 1 (which is the case 
of interest when the surface gravity k < 1), then Fig. shows clearly that 
-C 1 as well. Using small k approximations in ([13]) and ([19|) respectively 
and solving for k yields 

fctp « (Ulu) 1 ^ < 1, (20) 

consistent with our approximations. 

This very important result states that although the scale of the new 
physics is the lattice spacing 6(—l), the effects of the new physics occur long 
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before that scale is ever reached. With the ordinary wave equation the max- 
imum wavevector near the horizon is infinite due to the infinite blueshift 
(actually it is finite but trans-Planckian if the black hole is formed by col- 
lapse). One might have expected that on the lattice k tp would be of order the 
inverse lattice spacing 5 _1 but fl20D shows that this is not the case (although 
k ~ <5 _1 does occur far from the horizon — see for example the roots k + and 

in Fig. H] and the accompanying discussion). This fact — which is also 
true in continuum models with high frequency dispersion — was not noticed 
in earlier work on dispersive models. As long as uo <C 1 (|2~0D shows that the 
physics near the horizon that determines the Hawking flux depends only on 
the low order terms in k. This result is absolutely essential for the validity 
of the approximation used in the next section. 

In Section |6| we will discuss ways to avoid the problem of vanishing lattice 
spacing at infinity. This problem plays no role in the calculation of the rate 
of particle production however, so we will now explain how this rate can be 
obtained in a leading order approximation. 

5 Hawking radiation 

The lattice theory can be quantized in strict analogy with the quantization 
of linear field theory in curved spacetime so we will not spell it out here. A 
difference peculiar to the lattice theory (or dispersive continuum field theo- 
ries) is that the local notion of the ground state (or vacuum) is not Lorentz 
invariant but refers to the preferred free-fall frame. In a region where the 
function v(x) is constant — or is approximately constant on the scale of the 
relevant wavelengths — the line element (|3]) is flat and the action (|§D is that 
of a chain of identical masses coupled by identical springs. The ground state 
of this system is just the usual ground state of the normal modes, i.e., it is 
annihilated by annihilation operators for complex solutions to the oscillator 
equation with time dependence of the form exp(— iuigt) with positive c<j ff , that 
is, positive free-fall frequency (0). This is the free-fall vacuum. 

Given this initial vacuum state we would like to compute the particle flux 
seen by an observer sitting at a fixed location (fixed x coordinate) far outside 
the black hole. The natural notion of particle for such an observer coincides 
with that defined by Killing frequency, therefore we shall compute the number 
expectation value for an outgoing positive Killing frequency packet in a state 
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which at some initial time is the free-fall vacuum. The standard method 



of computing this|14| is to propagate the outgoing packet backward in time 
to the hypersurface where the vacuum state is defined. The norm of the 
negative free-fall frequency part of this packet is then (minus) the number 
expectation value. The norm referred to here is given by 

= iJ2 \v m (m<f>*m(t)d t (f) m (t) - Cj> m {t)d t( j>* m {t)), (21) 



and is the sum over a constant t surface of the t-component of the current 
associated with phase invariance of the action (||) (generalized to complex 
fields). 

Several methods can be used to compute the rate of Hawking radiation. 
One approach is to evolve a wavepacket backwards in time by numerical 
solution of the lattice wave equation (9), as was done by UnruhJTTJj. Al- 



ternatively, since the problem has time translation symmetry, one can just 
work with modes of definite Killing frequency. This is the approach we take 
here. The outgoing wavepacket is composed of wavevectors around k +s (and 
has positive Killing frequency) and arises from a pair of packets composed of 
wavevectors around k + and k_ respectively (which have positive and nega- 
tive free-fall frequency respectively). Using the arguments in 0, modified to 
the lattice model, it is straightforward to show that, for an outgoing packet 
narrowly peaked about the frequency u, the number expectation value is 

N(lA |(M^) + ^K(M^)AM 2 I (22) 

1 ) \(k +s (Lu)+Lu)v g (k +s (uj))c +s ^) 2 \ ^ ' 

where c_(c +s ) is the constant coefficient of the k_(k +s ) mode located far 
outside the black hole (where v (x) is essentially constant). We now turn to 
the computation off these coefficients. 



5.1 Mode equation 

The mode solutions to the lattice wave equation ([]) are of the form (|TT|), 
(0) 

<j) m (t)=e- wt f(m-t), (23) 

where to is the conserved Killing frequency. Plugging this into the equation 
of motion (||) produces a delay-differential equation (DDE), 

«(0 (/"(O + <2«/'(0 - ^ 2 /(0) + v'(0 (AO + 
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2 (/(g-i)-/(o) 2(/(p-/(e + i)) 

(T,(e-l)+i;(0) KO + «(e+l)) 







(24) 



where we have defined the new variable £ := (m — t), and := v(x(£)). A 
wavepacket that is outgoing at late times is composed of mode solutions that 
decay inside the horizon (see |§ for a discussion of the analogous boundary- 
condition in a dispersive continuum model). We therefore need to solve (p4|) 
subject to this boundary condition. 

The DDE ([2"4]) can be solved numerically, however it is more instructive, 
and sufficient for our purposes, to find an approximate analytic solution. We 
use the same analytical techniques as used in || . We first find an approximate 
solution (satisfying the above boundary condition) in a neighborhood of the 
horizon by the method of Laplace transforms, and then extend this solution 
far outside the black hole by matching to the WKB approximation. The 
mode coefficients q can then be read off directly. 

5.2 Near horizon approximation 

To solve the mode equation (p4j) near the horizon we first approximate v(£) 
as 



where k is the surface gravity of the black hole, and neglect all terms of order 
(k£) 2 . This requires that we stay close enough to the horizon that ^<1, 

Next we "localize" the DDE by first Taylor expanding /(£ — l),i>(£ — 
1), etc., and then truncating the expansions. Which terms to keep can be 
estimated as follows. The Taylor expansions produce the equation 



where we have grouped together terms in the expansion according to the total 
number of derivatives. The ellipses that appear inside parentheses denote 
other terms with a total of four derivatives and the other ellipses denote terms 
with six or more derivatives per term (only even numbers of derivatives occur 
in the expansion). Truncating the equation to second order in derivatives 
produces the ordinary wave equation. This is not sufficient for us because 



v(0 «-!+«£, 



(25) 



+ 



( 




+ ••• (26) 



14 



arbitrarily short wavelengths appear in the ordinary wave equation solution 
for the outgoing modes, so we must keep at least some of the higher derivative 
terms. 

Let us define an effective local wavevector fc(£) by /'(£)//(£) = 
Dropping the f^ m '(^) term compared to the f( iv '(£) term is accurate pro- 
vided that \k(£)\ <C 1 in the near horizon region |£| <C 1/k. We can es- 
timate from the dispersion relation in the near horizon approximation 
just as we did in section (f|). Outside the classical turning point (where 
2/3 /k), but still in a region where £ <C all relevant wavevectors 
are real and the largest wavevector behaves as ~ a/k£, and therefore sat- 
isfies \k(£)\ <C 1. For |£| < £ tp , the relevant wavevector becomes complex and 
has a magnitude \k(£)\ ~ cu 1 ^ 3 , therefore \k(£)\ <C 1 provided we only con- 
sider Killing frequencies satisfying to 1 / 3 <C 1. Even deeper inside the horizon 
where <C £ < — £ tp , the wavevector is approximately imaginary with 



magnitude again given by ~ y/c|£|, and therefere <C 1. Ignoring 

sixth and higher order derivatives in the equation (^) therefore requires that 

LU 1/3 < 1. 

To further simplify the equation, note that the ratio of the f^ v ' term to 
the /'" term is 



^ (27) 



/"'(£M0 

From above we know that|] fc(£) > c<// 3 , so we will have k <C provided 
that (jj 3> /t 3 . As long as is not ultra small therefore we need only keep the 
fourth order derivative termF] in the expansion (^). We therefore arrive at 
the ordinary differential equation (ODE) 

LfW _ 2<r _ 2(iw _ K)f _ _ K)f w (28) 

We show below by explicit calculation that the solution to (|28|) of interest to 
us is consistent with the approximations made above and therefore that this 
truncation of the mode equation is valid. 

9 Actually the wavevector of the outgoing wavepacket is smaller than this. For the 
outgoing packet though, all higher order derivative terms are negligible outside the classical 
turning point. 

10 We could in principle keep the third order derivative term as well and therefore enlarge 
the range of validity of our approximations in u>, however for simplicity we work with the 
simpler equation. 
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The ODE is the same as that considered in || (except for the coeffi- 
cient of the f^ v > term) where it was solved by the method of Laplace trans- 
forms with the same boundary conditions as discussed above. We therefore 
refer the reader to for the details of this computation. Using the saddle 
point approximation to evaluate the Laplace transform for £ 3> 1, we find 
that the solution satisfying the given boundary conditions can be expressed 
as 

/(0=/+(0 + /-(0 + /+.(O (29) 



where 



/+(0 « ^ e 3™/(2*)£-3/4-W(2«) exp fi|V24Kf 3/2 ) (30) 
/_(£) w jv e W(2«)^-3/4-W(2*) exp V24^£ 3/2 ) (31) 
/+,(£) » 2e m " /K sinh(7rcj//«)r(-zcj/fi;)f Lj/K (32) 



and 



N := e i7r / 4 ^(6K) 1 / 4 (V24^) ' . (33) 
To check the validity of our localization procedure, note for example that 

/+(0 /73 , .« 



/+(0 V V4 2^7 £ 



-W24< . (34) 



The absolute values of the two terms on the right-hand-side of (|34"D are both 
much less than one provided we restrict £ to the range 

K £ < k -1 (35) 

which was already assumed in making the saddle point approximation ( p5|) 
and the near horizon approximation (|25|) . Expression (|34]) is also in agree- 
ment with our earlier estimates of /'(0 obtained by estimating the position 
dependent wavevector Similar relations hold for the /_(£) and /+ s (£) 

modes as well. 



5.3 Match to the far zone 

The next step is to propagate the mode (|29|) away from the horizon to the 
constant v(£) region. This is accomplished by computing approximate solu- 
tions to the non-local DDE ( f24"|) by the WKB method. (Since the wavevectors 
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grow to order unity as v(£) goes to zero, we must use the full non-local DDE 
at this stage.) Some details of this computation are given in appendix |B| The 
result is that there exist three different WKB solutions which, when evalu- 
ated near (but not too near) the horizon, take the same functional forms as 
the Laplace transform solutions given by ([301, 31 ,[32 ). An appropriate linear 



combination of these WKB solutions can therefore be matched to the near 
horizon solution (p9l) yielding 



/(0 = v / 2^(e 3 " /(2K) /f KB (0 + e ™/«fWKB^ 

+ 2e™/«suih(™)r(-i^)fZ KB (0- (36) 



\ K J \ K 

Since the WKB approximation holds far outside the horizon, we are free 
to evaluate the solution there, and thus read off the constant coefficients of 
the modes exp(i/c£) with k = k +s , k + , fc_ in the constant v(x) region. These 
coefficients are simply given by the coefficients of the WKB solutions in (|36|) 
except that f^ KB also contain the amplitude factors (a;±27r) _1//2 respectively 
(see the appendix |B]). 

5.4 Kinematic factors 

The only remaining ingredient in evaluating the number expectation value 
(HH) is to compute the kinematic factors (k(u) + uj) and group velocity v g 
for each wavevector. From the dispersion relation (fl"3| ) (with the plus sign 
for the roots k +s , k + , k- corresponding to A, D and E respectively in Fig. |^) 
and the expression for the group velocity given by (|I3|) it is straightforward 
to show that 

(*H + u)v 9 {u) = ma(k ™~ IVI Mk/2). (37) 

Plugging in the small \v\ expressions for the k_ and k +s wavevectors com- 
puted in Section [|, we find that Q3"7| ) reduces to —{uj — 2h)/\v\ for the k_ root 
and oj/\v\ for the k +s root. Putting all these results together we find for the 
number expectation value (P2|) 



17 



where T# = k/2tt is the Hawking temperature. Therefore we see that to 
leading order in the lattice spacing the particle flux is thermal at the Hawking 
temperature in agreement with the ordinary wave equation. 

This derivation is valid as long as (i) k <C uj 1 ^ <C 1 and (ii) the WKB 
approximation can be used to connect the far zone with the zone <C 1 
near the horizon. This last condition should be satisfied as long as to is 
not extremely small compared to k, although we shall not attempt to write 
out the general conditions here (which are possibly more restrictive than the 
k <C to 1 / 3 condition already given). 

6 Models with finite lattice spacing at spatial 
infinity 

One way to avoid the problem of vanishing lattice spacing at infinity is to 
simply not let v(x) go to zero at infinity. It might seem that we have no 
freedom to make this choice, since the asymptotic form of the metric is 
determined by the black hole. However, we need not use a free-fall coordinate 
that is at rest at infinity. Instead, the coordinate lines can be chosen moving 
uniformly toward the black hole at infinity. In appendix |A] it is shown that, 
in terms of the proper time t' along the congruence of infalling geodesies of 
energy E > 1 and the proper distance x' along the spacelike slices orthogonal 
to these geodesies, the line element takes the form dt' 2 — (dx' — VE(x')dt') 2 
for some function ve- Note that this is the same form as ([I]), with a different 
function uj^u which, in particular, does not vanish at infinity: Ve{oo) = 
— [E 2 — l) 1 / 2 . Proceeding as before one then arrives at the new line element 
(|3p, but with v replaced by ve- With this choice the preferred frame is not 
asymptotically at rest with respect to the black hole. Although this certainly 
solves the problem from a mathematical point of view, it is not physically 
satisfactory. Our "in" vacuum boundary condition depends on the choice 
of the preferred frame, and it just doesn't make much sense to rely on the 
assumption that the black hole is moving relative to the vacuum. 

A more satisfactory resolution would be to choose the discretization such 
that the lattice spacing is a fixed proper distance on some initial slice. If we 
then let the lattice points fall into the black hole, the proper lattice spacing 
will not remain constant on the surfaces of equal proper time. Nevertheless, 
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such a lattice will be perfectly well behaved at infinity, and the time depen- 
dence will be invisible to long wavelength modes that do not "see" the lattice 
at all. Although they are not ultimately satisfactory, we think it is instruc- 
tive to understand the physics of such models with growing lattice spacing. 
We now describe a class of such models. 

It is only necessary to reparametrize the y coordinate ([J) before discretiz- 
ing. To this end, we define a new coordinate z by 

W{z) = y = t- [ X dx'/v(x'), (39) 

where Xh is the value of x at the event horizon, i.e., v (xh) = — 1- The original 
x coordinate measures proper length on a constant t surface in the metric 
([!]), so we choose z to agree with x at t — 0. This implies 

W(z) := - f dx'/v(x'). (40) 

JXh 

In terms of the function W, the defining relation for z can be written as 

W(z)=t + W(x), (41) 
which can be solved for x(t, z) as 

x = W~ l (W(z) - t) . (42) 
In the coordinates (t, z) the line element (HD becomes 

ds 2 = dt 2 - ( V -^l\ 2 dz\ (43) 
\v(z)J 



where x(t, z) is the function defined by (f42j). In these coordinates the Killing 
vector x (which is d t in the (t, x) coordinates and d t + d y in the (t, y) coor- 
dinates (^)) is given by 

X = d t - v(z)d z . (44) 

When d z is modified in the action, either by higher derivatives or discretiza- 
tion, the presence of the factor v (z) in ([44]) will prevent the survival of the 
symmetry generated by £. Not even a discrete remnant of the symmetry 
survives in the discrete case. 
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At any finite t, the the spatial scale factor v(x)/v(z) goes to unity as z 
goes to infinity, as long as v(x) goes to a constant (including zero) at infinity. 
Thus, the coordinate z always measures proper distance sufficiently far from 
the black hole. Along a line of fixed z, v(x)/v(z) grows as a function of t 
as the horizon is approached, since x is getting smaller and we are assuming 
\v(x)\ grows as x decreases. That is, the proper spacing of the z coordinate 
grows with t because of the relative acceleration of the free-fall worldlines. 

At the horizon v(xh) — — 1, W(xh) = 0, and therefore z = W / ~ 1 (t). This 
yields the form of the line element evaluated at the horizon: 



7 2 1 i 2 /.,■_. i . . \ 2 

horizon ^ 



(w-\t))] 2 dz 2 (45) 



Let us now consider two examples to see what this coordinate change 
yields. First, consider the Schwarzschild line element, for which v(x) = 
— (2kx) -1 / 2 , where k is the surface gravity 1/AGM. In this case the line 
element ([43D becomes 

^-^-{'-liSwT' 3 ^ < 46) 

and at the horizon this reduces to 

ds 2 |horizon = dt 2 - (1 + 3nt) 2/3 dz 2 . (47) 

For numerical calculation, it would be more convenient to have v(x) go to 
zero more quickly than x -1 / 2 , so let us also consider the exponential velocity 
v (x) — — exp(— kx). In this case the line element ([4"3|) becomes 

ds 2 = dt 2 - (l - Kte~ KZ ) ~ 2 dz 2 , (48) 
and at the horizon this reduces to 

ds 2 (horizon = dt? ~ (1 + Ktf dz 2 . (49) 

Discretizing the z coordinate will yield a new lattice theory in which the 
proper lattice spacing is constant at infinity, so it is possible to propagate 
wavepackets in a sensible way all the way out to where v(z) m 0. Therefore 
the ingoing waves that produce the outgoing waves must originate at infinity 
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as combinations of the standard flat space lattice modes. No exotic low 
frequency modes are available in this case. The low frequency ingoing waves 
behave like ordinary continuum ingoing waves which sail right across the 
horizon. They will not bounce off the horizon. So where can a low frequency 
outgoing mode come from? 

The lack of even a discrete time translation symmetry seems to provide 
the answer. When a low frequency outgoing wavepacket is propagated back 
close to the horizon, it gets blueshifted. Eventually its wavevector gets so 
large that it can sense the lack of time translation invariance in the lattice the- 
ory. At that point, there is no longer any reason for its Killing frequency to be 
conserved. Using an eikonal approximation we will show in the next section 
that the Killing frequency is indeed shifted so that, when the wavepacket 
propagates backwards in time back out to infinity, it arrives with a large 
wavevector, on the order of the lattice spacing, and a correspondingly large 
Killing frequency. At this stage we have no solid proof that waves on the 
z-lattice will behave in the way indicated by the eikonal approximation. It 
should be possible to adapt Unruh's numerical computation on the y-lattice 
to see what in fact happens on the z-lattice. 



7 Origin of the outgoing modes revisited 

In deriving the eikonal approximation we forget that the space is discrete 
and just make the substitution d z — > exp(<9 z ) — 1 in the continuum action in 
(t,z) coordinates (in units of the lattice spacing): 




This leads to an infinite order PDE to which the standard eikonal or geomet- 
rical optics approximation can be applied. One assumes that the wavelength 
and period of the wave are short compared with the length and time scales 
on which the background is varying and slowly changing on their own scales. 
This is reasonable for much of the trajectory of the wavepackets we are inter- 
ested in, but the latter condition fails at the turning point near the horizon. 
Nevertheless, the results obtained in this way seem reasonable and we would 
be surprised if a lattice calculation failed to confirm the general picture pro- 
vided by this approximation. 



21 



Making this approximation, and assuming a wave of the form 

exp(— iuit) exp(ikz), (51) 
we arrive at the dispersion relation 

u 2 = -g zz {t,z){F{k))\ (52) 
where the function F(k) is given by 

F[k) = 2 sin(/c/2) (53) 

and, using (pf), 

-g*%z) = [v(z)/v{x(t,z))] 2 . (54) 

Note that now uj (rather than uig) stands for the free-fall frequency. 

The eikonal approximation in this case amounts to Hamilton's equations 
for the phase space variables (z, k) with the Hamiltonian 

H = ± v ^g~^F(k). (55) 

The free-fall frequency is just uj = H, so the sign of this frequency is the 
sign of ±F(k). On the lattice, wavevectors differing by 27m are identified, 
so a complete set of k values is the interval [0, 27r). In this interval F(k) in 
( |53|) is positive, so the sign of uj is the sign of the prefactor ±. Instead of 
keeping this prefactor alternative, we can double the range of k to (— 2tt, 2tt) 
and always use the + sign in the Hamiltonian (pH|), since F(k) = —F(—k) is 
negative when k G (— 2ir, 0). 
Hamilton's equations are 

dz/dt = d k F (56) 

dk/dt = -d zS [=g^F. (57) 

We have solved these equations numerically for the case of the exponential 
velocity function v(x) = — exp(— kx), for which (|54]) yields 

^J-g zz (t,z) = 1 - Klzr**. (58) 

We used k = 0.001 and started the trajectories at the initial position z(0) = 
10, 000 at t — 0. The unit here is the lattice spacing in the ^-coordinate. For 
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each initial wavevector k(0) we obtain a trajectory (k(t), z(t)). To visualize 
the results, it is convenient to plot k(t) versus v (x(t)) because the value of 
v(x) indicates the static radial position whereas the z coordinate lines are 
falling. (We could also have plotted versus x(t) itself but it is helpful to be 
able to see the value of v(x) on the same graph.) The results are given in 
Fig. |3|. The equations of motion fl56|) and (|57|) are symmetric under k — > —k, 
so the solutions for negative k's are obtained by changing the sign of k. 



6.0 



k 




Figure 3: Plot of the wavevector trajectories as a function of the background 
free-fall velocity function v (x). 

At spatial infinity, where v (x) = 0, the right moving modes have k 6 (0, 7r) 
and the left moving modes have k e (tt,2tt). Thus we send in modes with 
k in ("7T, 2tt). The ones near 2tc are equivalent to ordinary small negative k 
modes and just cross the horizon. Since the group velocity (j56|) is always less 
than or equal to the speed of light (—g zz (dz/dt) 2 = (e\F) 2 = cos 2 (A;/2) < 1), 
these modes can never return to the outside once having crossed the horizon. 
Coming down from 2tc, at some critical value of k there is a trajectory that 
asymptotes to the horizon and zero k. Below this critical k are the exotic 
modes that bounce off the horizon and return to spatial infinity. The crucial 
thing to notice here is that an exotic ingoing mode can produce a non-exotic, 
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very low wavevector outgoing mode. This is only possible because the lattice 
equations violate time translation symmetry at short wavelengths, so there 
is no conservation of Killing frequency to prevent this from happening. 

8 Discussion 

It is intriguing that violation of time-translation invariance visible only at 
short wavelengths plays a crucial role in accounting for the outgoing modes. 
In our model this time-dependence is a consequence of the growing lattice 
spacing due to spreading of free-fall trajectories. At a more fundamental 
level, one expects the Killing symmetry of a black hole background to be vi- 
olated by the gravitational back-reaction to the quantum fluctuations of the 
matter fields. A vague suggestion was made in JT5] that the back-reaction 
might evade the conservation of Killing frequency and allow the outgoing 
modes to originate as ingoing modes from spatial infinity. Our simple model 
studied here seems to lend credence to this hypothesis, although the im- 
plementation is still in a background field approximation and has nothing 
obvious to do with the back-reaction. 

It is scary to be violating time-translation invariance in the lattice theory. 
However, the characteristic timescale is long, ft -1 according to either (|4"?| ) or 
([49]) for example, and even this time dependence is invisible to wavelengths 
long compared with the lattice spacing. It therefore seems that the low energy 
physics is immune from direct effects of this violation of time-translation 
symmetry, even though the outgoing modes owe their very existence to this 
violation! 

We still do not have a satisfactory discretization of field theory in a black 
hole background. Either our lattice spacing goes to zero at infinity, or it 
grows as points fall in towards the horizon. For the Schwarzschild metric, 
the total amount of growth during the Hawking lifetime M 3 is, from (|4"T|), 
of order M 2//3 . Thus if the lattice starts out with Planck spacing, it ends 
up with spacing of one Angstrom after the evaporation of a solar mass black 
hole. But this is only the radial spacing. If the lattice points are falling 
on radial trajectories from radius r-i to t\ their transverse proper spacing 
decreases by the factor ri/r 2 . 

It seems that to maintain a uniform lattice spacing in some preferred 
frame with a freely falling lattice of fixed topology is not possible. This 
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suggests that one should be thinking about a lattice in which points can be 
created or annihilated in order to keep the spacing uniform. 

An expanding cosmology provides a simpler setting than the black hole 
in which to contemplate the lattice question. As the universe expands, the 
lattice spacing will grow if the lattice points are at rest in the cosmic rest 
frame. Weiss |16| confronted this issue in trying to formulate lattice field 
theory in an expanding universe. He noted a very interesting point: if the 
couplings of an interacting field theory are fixed on the expanding lattice, 
then the renormalized parameters at a fixed proper scale will depend strongly 
on the cosmological epoch. One could of course adjust the lattice parameters 
as the scale factor evolves, but from a fundamental point of view that is 
artificial. Moreover, if the lattice spacing started out in the early universe at 
the Planck scale, it would quickly become too large to appear continuous at 
large scales. Both these problems would be eliminated if the lattice were itself 
dynamical, with points being added at the right rate to keep their density 
constant. 

Allowing the lattice topology to be dynamical thus seems very natural. 
It would be interesting to see if field theory can be sensibly formulated on 
dynamical lattice models and, if so, to study the consequences for cosmology 
and black hole physics. 
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A Free-fall coordinates 

In this appendix we show that in a general static two-dimensional spacetime 
coordinates can always be chosen (at least locally) so the line element takes 
the form ([!]). Let x a De the time-translation Killing field, let u a be the unit 
tangent vector to a congruence of timelike geodesies all of the same energy E 
and invariant under the symmetry, and let s a be the (unique up to sign) unit 
vector orthogonal to u a . Then u a = E~ x \ a + vs a , where v 2 = 1 — (x 2 /E 2 ). 



25 



The assumed symmetry of u a implies [x, s] = 0, so there exist coordinates r 
and x such that E~ l x a = {d T ) a and s a = (d x ) a . In these coordinates the line 
element takes the form 

ds 2 = (1 - v 2 )dr 2 + 2v drdx - dx 2 (59) 
= dr 2 - (dx - v dr) 2 . (60) 

Note that r coincides with the proper time along the orbits of u a , the lines of 
constant r are orthogonal to these orbits, and x measures the proper distance 
along these lines. Note also that, because of the symmetry, v(t,x) = v(x) 
depends only on the coordinate x. If \ is normalized at infinity we have 
v(oo) = ±(l-E- 2 ) 1 ' 2 . 



B WKB solutions to the DDE 

In this appendix we discuss the application of the WKB approximation to 
finding approximate solutions to the DDE (p4|) . We assume a solution of the 
form 

/(0=expf+i fm) (61) 



and substitute into the DDE This results in the equation 

v(0(+^(0 - k 2 (0 - 2uk(0 - u 2 ) + tv'(0(H0 + w) (62) 
2(exp(-z 1^ k(u)) - 1) + 2(l-expH// +1 fc( M ))) = Q 



We can rewrite the exponentials in a form more appropriate for the WKB ap- 
proximation by Taylor expanding the integrand about £ and then evaluating 
the integrals, eg., 

£ +1 k(u)du = k{i)+ l -k> {£) + ■■■. (64) 

For bookkeeping purposes, it is now convenient to make the substitution 
£ — ► a£, which has the effect of scaling n th order derivatives in the equation 
by l/a n . Now expand fc(£) as 

fc(0 = *b(0 + -*i(0 + -"> (65) 
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substitute into (|63|), and demand that each coefficient of the separate powers 
of 1/a vanish. The leading order equations are 



i' 2 (0(*o({) = (2sin(A„(fl/2))' (66) 

= in -5^). (67) 



The first of these equations is of course the dispersion relation flTBD that we 
derived in Section while the second produces the first order correction to 
the leading order root from the dispersion relation. 

To solve the dispersion relation near the horizon (where v ~ — 1) note 
that when u <C 1 then 2sin(A;o/2) « (k — /cq/24). Using this approximation 
it is straightforward to show that the roots are 

k 0>± (0 » ±v/24(l- KOI)- ^!^,) (68) 

feo 1+s - I3R ^. (69) 

Substituting these into the expression for k\ above gives the first order cor- 
rection term. 

To match the WKB solutions given here to the Laplace transform solu- 
tions given in Section [| we need only substitute the near horizon expansion 
for v(£) ~ — 1 + k£ into the expressions for ko and k\ and evaluate the inte- 
grals given in fl6~l|). Note that k\ will yield in general a non-trivial amplitude 
factor. 
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